/ 



} A 


NASA Tfechnical Memor^dum 1061^6^ 


h:. 




ProbabiliMic Simulat ion of .Malti-Scale 
jj^mposite BehayiQr ^ ^ ^ 










” D.G. tiaw, Nf.C. Shiao, and S.N. Singhal 
Sverdrup^echnology, Inc. _ 

Lewis Research Center Groi^ 

Brook Park, Ohio ~ 


and. 


Christo s C, C hamis _ _ . ; ' , 

NatioTial Aeronautics and Space Administration 


Cleveland, Ohio 






^ r Tg w g K aP i ^''.r:i' ; , .^4 i"'*. ■ ^ tv:' 




Prepared for the ■ • ■ - 

Sixth Japan-United States Conference on Composite Materials 

Orlando, Florida, June 22-25, 1992 : ; 


IM/V5A 


-- xmaca TM*-10fil96) probabilistic 

siwtlnol Of MULTI-SCALE COMPOSITc 

-i; behavior (NASA) lA p 


N93-28633 
Unci as 


J 


dm 


03/39 0170793 




PROBABILISTIC SIMULATION OF MULTI-SCALE COMPOSITE BEHAVIOR 


D,G. Liaw, M.C. Shiao, and S.N. Singhal 
Sverdrup Technology, Inc. 

Lewis Research Center Group 
Brook Park, Ohio 44142 

and 


Christos C. Chamis 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

A methodology is developed to computationally assess the probabilistic composite material 
properties at all composite scale levels due to the uncertainties in the constituent (fiber and 
matrix) properties and in the fabrication process variables. The methodology is 
computationally efficient for simulating the probability distributions of material properties. 
The sensitivity of the probabilistic composite material property to each random variable is 
determined. This information can be used to reduce undesirable uncertainties in material 
properties at the macro scale of the composite by reducing the uncertainties in the most 
influential random variables at the micro scale. This methodology has been implemented 
into the computer code PICAN (Probabilistic Integrated Composite ANalyzer). The accuracy 
and efficiency of this methodology are demonstrated by simulating the uncertainties in the 
material properties of a typical laminate and comparing the results with the Monte Carlo 
simulation method The experimental data of composite material properties at all scales fall 
within the scatters predicted by PICAN. 



INTRODUCTION 


Polymer matrix composites (PMC) form an important class of engineering materials 
in structural applications. Their outstanding mechanical pn ^)er t ies , such as durability and 
corrosion-resistance characteristics, combined with low density are very attractive to the 
aerospace industry. The mechanical properties of composite materials depend on a wide 
variety of variables, for example: the fiber and matrix material properties; and fabrication 
variables at all stages of fabrication process, such as, fiber volume ratio, void volume ratio, 
fiber misalignment and ply thickness. These variables are statistical in nature, therefore, 
the mechanical properties of a composite material should be quantified probabilistically. As 
a consequence, the behavior of composite structures shows a scatter from its average value. 
Traditionally, an ad-hoc "safety factor" is used in the design to account for the difficulty 
in predicting the structural behavior. However, this approach may result in either an 
ultraconservative or an unsafe design. 

In order to probabilistically assess the behavior of composite structures, a 
methodology was developed at the NASA Lewis Research Center for simulating the 
probabilistic material properties in all composite scales [1-3]. The methodology requires 
the idCTtification of uncertainties at constituent and ply levels including material properties 
and fabrication variables. The random variables are then filtered through an analyzer which 
combines the composite mechanics and probability theory to computationally simulate the 
probabilistic material properties. Although the uncertainties in material properties can be 
computed by the Monte Carlo simulation (MCS) method, it is inefficient and expensive. 
In order to save computational time, the newly developed methodology uses the Fast 
Probability Integration (FPI) algorithm, instead of the traditional Monte Carlo simulation, 
to provide an efficient and affordable way for computationally simulating the uncertainties 
in composite material properties. The methodology was implemented in the computer code 
PICAN (Probabilistic Integrated Composite ANalyzer) (Ref. [4]) and is described in the 
present paper. The uncertainties of composite material properties at ply and laminate levels 
were simulated via PICAN. The sensitivity of a given material property to each random 
variable was also quantified. 


PROBABILISTIC SIMULATION OF COMPOSITE MATERIALS 

Composite materials are made of fiber and matrix through an appropriate fabrication 
process. Since the mechanical properties of the constituent material vary statistically, the 
uncertainties in composite material properties are produced from the uncertainties in the 
constituent material properties, also from the uncertainties in the fabrication parameters. 
The probabilistic simulation of composite material properties requires identification of 
uncertainties at all composite scales as described below. 

Many thermoset laminated fiber-reinforced composite materials are manufactured by 
three basic steps: tape, layup and curing (see Fig. 1). The uncertainties incurred during the 
layup process are due to the misalignment of ply orientation. Typical uncertainties incurred 
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from the curing process are intralaminate voids, incomplete curing of resin, excess resin 
between plies, excess matrix voids and porosity, and variation in ply thickness. 
Uncertainties incurred in the constituent material properties include all the fiber and matrix 
material properties. These uncertainties are referred to as the random variables (as 
enumerated in Figure 1) which will be processed computationally through composite 
mechanics and probability methods as depicted in Figure 2. The details of the computer 
code PICAN for the computational simulation of uncertainties in composite material 
properties are described next. 


PROBABILISTIC INTEGRATED COMPOSITE ANALYZER ; PICAN 

PICAN (Probabilistic Integrated Composite ANalyzer) is a computer code for the 
probabilistic simulation of composite material properties. It integrates two NASA in-house 
computer programs, ICAN (Integrated Composite ANalyzra) (Ref. [5]) and FPI (Fast 
Probability Integrator) (Ref. [6]). ICAN was developed based upon composite mechanics 
to simulate the composite material behavior at all composite scales including many effects 
of environment and fabrication. FPI, developed using structural reliability algorithms, 
probabilistically processes the material property information generated by ICAN and 
efficiently computes the cumulative distribution functions (cdf) of composite material 
properties at all composite scales. Further details of ICAN and FPI are explained in the 
following. 

ICAN MODULE 

The ICAN module [5] performs a comprehensive analysis of composite material 
properties at all scales by selectively perturbing the random variables at the constituent level 
(fiber and matrix material properties, fiber volume ratio and void volume ratio), and the ply 
level (ply misalignment and ply thickness). Perturbed material properties at all composite 
scales corresponding to all ^e perturbed random variables are computed and stored in a 
database. 

FPI MODULE 

The FPI module [6] performs the probability calculations. It takes the relationship 
between the composite material property and the random variables as well as the probability 
distributions of the random variables to compute the probability distribution of that material 
property. The relationship is expressed as 

1-1 i -1 


where Mp is any material property; N is the number of random variables considered; Vj are 
the random variables; ao, ai, bj are calculated constants. More theoretical details concerning 
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FPI can be found in Ref. [6]. 

PICAN SIMULATION 

The random variables at the constituent and the ply levels are first identified and 
their respective distribution types and corresponding parameters are predetermined. These 
variables are then selectively perturbed in order to create a da t a ba se to determine the 
constants a„, aj, and bj in the approximate relationship giv«i by Eq. (1) between the 
composite material property and the random variables. For each set of perturbed values of 
random variables, composite micro-mechanics is used to determine the corresponding 
mechanical properties at the ply level. Laminate theory is thrai used to determine the 
corresponding mechanical properties at the laminate level. The process is repeated until 
enough data sets are generated by ICAN module that the functional relationship between 
material properties and random variables can be numerically determined through a least 
square curve fitting procedure in FPI module. For example, if each random variable is 
perturbed four times, a quadratic functional relationship can be obtained. With the known 
probabilistic distributions of the random variables and the numerically determined 
relationships between the composite material properties and the random variables, the fast 
probability integration is applied. For every discrete response value, a corresponding 
cumulative probabiUty is computed by FPI. This process is rq)eated un^ the cumulative 
distribution function is represented with sufficient accuracy. The output information from 
FPI for a given material property includes its discretized cdf values, the coefficients for a 
special type of probability distribution function, and the sensitivities of the composite 
material property to the random variables. 


NUMERICAL RESULTS AND DISCUSSIONS 

The methodology described previously was applied to a 16-ply (O 2 / ±45/02/90/0)s 
graphite/epoxy (AS/3501-5) composite with 1.8% moisture at room temperature. The 
experimental data for the longitudinal modulus, transverse modulus, in-plane shear modulus, 
major Poisson’s ratio, longitudinal and transverse tensile strengths, longitudinal and 
transverse compressive strengths, and in-plane shear strength of the graphite/e^xy 
unidirectional single ply and 16-ply composite laminate can be found in Ref. [7] and listed 
in Table 1. Based upon these experimental results (mean values) for the ply and laminate 
properties, the mean values of constituent material properties (fiber and matm) were 
estimated. The assumed statistical properti^ (mean value, coefficient of variation (coy) 
or standard deviation (stdv), and distribution type) of each random variable are listed in 
Table 1. Typical probabilistic material properties for ply and lairunate moduli and strengths 
were simulated. The sensitivities of the ply and laminate properties to the most sensitive 
random variables were also computed at two cumulative probability levels, 0.(X)1 and 0.999. 
All the cumulative distribution functions of the material properties at ply and laminate levels 
were compared with the results predicted by a Monte Carlo simulation with 1000 samples. 
Good agreement was observed. 
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The cumulative distribution functions (cdf) of the ply moduli and the sensitivity 
information are shown in Figures 3 to 5. In each Figure, the experimental datum (mean 
value) is indicated by an arrow. In Figure 3, it is fotmd that the ply modulus in the fiber 
direction (E,u) is most sensitive to the fiber modulus and fiber volume ratio. This means 
that the reduction of the uncertainty in the ply modulus can be achieved mostly by the 
reduction of the uncertainties in the fiber modulus and fiber volume ratio. In Figure 4, it 
is found that the shear modulus of matrix and fiber volume ratio have the most contribution 
in the uncertainty of the ply in-plane shear modulus (£,, 2 ). The uncertainty in the ply major 
Poisson’s ratio (i/,i 2 ) results mainly from the uncertainties in the fiber and matrix Poisson’s 
ratios and the fiber volume ratio as shown in Figure 5. The probabilistic tensile strength 
in the fiber direction (S,i,x) was simulated and shown in Figure 6. From the simulation, the 
probabilistic tensile stren^ is most sensitive to the uncertainties in the fiber tensile strength 
and fiber volume ratio. 

The cumulative distribution functions (cdf) of the equivalent laminate moduli and 
their sensitivities to the random variables are shown in Figures 7 to 9. It is found that both 
the modulus in the X direction (Et^J and the in-plane shear modulus (Etxy) are sensitive to 
the uncertain fiber modulus and fiber volume ratio as shown in Figures 7 and 8. In Figure 
9, the uncertainty in the equivalent major Poisson’s ratio results from the uncertainties 
in fiber modulus, fiber volume ratio, fiber and matrix Poisson ratios, matrix shear modulus 
and ply misalignment. Based on Ae strengths calculated in the ply level, the laminate 
tensile strength in the X direction (Sc„t) was predicted by assumed both first ply failure and 
fiber fracture failure modes as shown in Figure 10. It is shown that neither scatter predicted 
by these two criteria covers the experimental scatter. However, the experimental scatter 
Mis between the scatter domain predicted by both failure criteria. Indeed, the assumed 
failure modes represent extreme "failure" sceneries in a composite specimen. Consequendy, 
the predicted cdf s describe upper and lower bounds of the laminate strength. From the 
sensitivity analysis, the tensile strength predicted by first ply failure mode is most sensitive 
to the fiber modulus and matrix tensile strength. On the other hand, the tensile strength 
predicted by fiber fracture mode is most sensitive to fiber tensile strength and fiber volume 
ratio. 


CONCLUDING REMARKS 

A methodology has been developed to accurately and efficiendy simulate the 
uncertainties in composite material properties at all scales of the composite considering the 
uncertainties in all variables describing the composite material behavior (constituent material 
properties and fabrication variables). This methodology, integrating the probabilistic 
concepts in conjunction with composite mechanics, has been implemented into the computer 
code PICAN. 

The accuracy and efficiency of the present methodology were demonstrated by 
comparison with the Monte Carlo simulation method. The uncertainties in ply and laminate 
properties were quantified. The experimental data of the composite material properties at 
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aU scales feU within the scattere predicted by PICAN. The sensitivity of the mateii^ 
property to each random variable at the constituent and ply levels was computed. This 
information can be used as a guide to increase the structural reliability or to reduce the cost. 
Furthermore, the probabilistic structural analysis and risk assessment can be pCTformed once 
the uncertain ply and laminate prt^jerties are computationally simulated. It is also 
interesting to note that the most sensitive random variables for the laminate strengtiis 
predicted by different failure mode may be different. 
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TABLE 1- STATISTICAL PROPERTIES OF RANDOM VARIABLES 


Random Variables 

Ref. [4] 

Mean 

Value 

cov (or 
stdv) 

Distribution 

Type 

Rber: 





Normal modulus Em 

32 mpsi 

32 mpsi 

8 % 

Normal 

Normal modulus £^22 

2 mpsi 

3 mpsi 

8 % 

Normal 

Poisson’s ratio V 12 

0.2 

0.23 

8 % 

Normal 

Poisson’s ratio ^23 

0.2 

0.25 

8 % 

Normal 

Shear modulus Gh 2 

2 mpsi 

2.5 mpsi 

8 % 

Normal 

Shear modulus G {23 

2 mpsi 

2.5 mpsi 

8 % 

Normal 

Tensile strength Sn- 

400 ksi 

400 ksi 

8 % 

Weibull 

Compressive strength Sr; 

400 ksi 

400 ksi 

8 % 

Weibull 

Matrix: 





Normal modulus 

0.4 mpsi 

0.45 mpsi 

8 % 

Normal 

Poisson’s ratio 

0.4 

0.41 

8 % 

Normal 

Tensile strength S„,x 

7 ksi 

6.7 ksi 

8 % 

Weibull 

Compressive strength S„^ 

36.3 ksi 

39 ksi 

8 % 

WeibuU 

Shear strength S^ 

7 ksi 

8.9 ksi 

8 % 

Weibull 

Fabrication Variables: 





Fiber volume ratio (fvr): 

60% 

60% 

8 % 

Normal 

Void volume ratio (wr): 

0 . 01 % 

0 . 01 % 

8 % 

Normal 

Ply thickness: 

0.0055 in. 

0.0055 in. 

5% 

Normal 

Ply misalignment: 

0 

0 

0.9® stdv 

Normal 
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Figure 1. Random Variables in the Typical Fabrication Process 
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Figure 2. Probabilistic Simulation of Composite Materials 
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e 3. CDF’s and Sensitivity Factors of Ply Modulus in X-Direction 



Figure 4. CDF’s and Sensitivity Factors of Ply In-Plane Shear Modulus 
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Figure 6. CDF’s and Sensitivity Factors of Ply Tensile Strength in X-Direction 
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7. CDF’s and Sensitivity Factors of Laminate Modulus in X-Direction 



Figure 8. CDF’s and Sensitivity Factors of Laminate In-Plane Shear Modulus 
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MAJOR POISSON’S RATIO SENSITIVITY FACTORS 

Figure 9. CDF’s and Sensitivity Factors of Laminate Major Poisson’s Ratio 
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